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1. Introduction 

The possible equivalence of string and gauge theory has long history going back to birth of 
string theory as a possible theory of the strong interactions and the recognition that Feynman dia- 
grams of perturbative large N gauge theory could be naturally classified into a sum over surfaces 
of given genus. 

However, there has been renewed interest in the subject in recent years stemming from Malde- 
cena's discovery that the type II string theory in anti-deSitter space could be described in terms 
of a particular supersymmetric gauge theory in four dimensions [[]]]. This theory has sixteen real 
supersymmetries and gauge group SU (N). This so-called AdSCFT correspondence has been tested 
most extensively in the limit N — > °° and for low energies in which case the string theory reduces to 
classical supergravity. However it is widely conjectured to be true away from these limits. Perhaps 
surprisingly the gauge theory in question looks very different from QCD - it is a superconformal 
field theory with neither asymptotic freedom nor confinement. 

Since then many other such dualities have been postulated corresponding to varying the di- 
mension of the gauge theory, adding additional matter fields and breaking all or part of the su- 
persymmetry and conformal invariance. The dual gauge systems possess the common feature of 
residing on the boundary of the space on which the string theory is defined and are hence often 
referred to as giving a holographic representation of the gravitational system. One such class of 
dualities arise naturally from the original AdSCFT construction and consist of a mapping between 
type II string theory containing N Dp-bmnes and (p+l) dimensional supersymmetric gauge theo- 
ries with gauge group SU (N) [Q]. In this talk we will concentrate on perhaps the simplest of these 
systems which describes the dynamics of DO-branes in string theory in terms of super Yang Mills 
quantum mechanics [|[]. 

We will be interested in this mapping at finite temperature corresponding to a background 
string geometry in which the time coordinate is Euclidean and periodic. At low energies the resul- 
tant supergravity equations contain a black hole with charge N and metric 



ds 1 = a'[-h{U)dt z + h-\U)dU 1 + — —d&l (1.1) 

U 2 



where (A = Ng s ot' 3 ^ 2 ) and the function h(U) is given by 



h(U) = -^7= ^-[^) 1 <1.2) 




C2 VX 

with a' = ^ tring is the (inverse) string tension and g s the string coupling. Such a black hole has a 
temperature and entropy given by 

The holographic conjecture states that the dual Yang-Mills model has ,/V colors and is to be com- 
puted at the same temperature with X identified as usual 't Hooft coupling A = g\ M N. The state- 
ment of duality implies that the free energies of both gauge theory and black hole are equal. Duality 
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thus offers a way of understanding black hole entropy as arising from a counting of microstates in 
the dual gauge theory. 

The analysis that leads to the duality conjecture is only valid in which the string theory reduces 
to supergravity. This requires Ng s » I and a' — ► 0. Thus the dual super Yang Mills quantum 
mechanics is strongly coupled. This is rather a generic feature of these gauge-gravity dualities - 
typically the regime in which the string theory is tractable corresponds to a strongly coupled large 

./V gauge theory. Furthermore, in the quantum mechanics case the dynamics of the gauge theory 

i 

depends only on the dimensionless coupling j8 = ^ which then implies it is the low temperature 
behavior of the gauge theory that provides a description of the dual black hole thermodynamics. 

It is not known what happens as we increase the temperature of the gauge theory. Familiarity 
with thermal gauge theories leads us to speculate that the system can undergo a deconfining phase 
transition. Indeed, as we shall show just such a transition occurs in the quenched theory. In the 
gravitational language raising the temperature corresponds to increasing a' and with it the strength 
of classical string theory corrections. It is not known for sure what happens in this case. It is 
possible that the duality breaks down and no correspondence exists between the gauge and gravity 
models. Or more exotic possibilities could exist with the existence of a thermal phase transition in 
the gauge theory signaling some sort of phase transition in the gravitational system - for example 
a transition from black holes to a hot gas of strings and branes. One of the goals of our work will 
be to map out this gauge-gravity correspondence and to understand the phase diagram of the string 
theory by conducting numerical simulations of the gauge theory. 

In this talk we will report on numerical simulations of a related model with four supersym- 
metries. General arguments exist [Q] that this model may lie in a similar universality class to its 
sixteen supercharge cousin. This theory has also been studied recently in 



2. Action and supersymmetries 

The £? = 16 super Yang-Mills QM is obtained by dimensional reduction of jV = 1 SYM in 
D = 10 dimensions. 

S = j J R dT ((D T X ; -) 2 - [X^Xj] 2 + 2m> a D^ a +2¥°(T i )5[X«,'Fp]) (2.1) 

where x V a is a 16-component spinor and T, = jojj where i = . . . 10 are the D = 10 Dirac matrices 
in Majorana representation. The bosonic sector of the model consists of 9 scalars Xj and a gauge 
field A. All fields take values in the adjoint representation of the gauge group. The supersymmetries 
are given by 

8A = -2i x ¥ a e a (2.2) 
5Xi = -2e"(r i )^ / 3 (2.3) 

8v a = \ (F,)fiD t x t + i([r ; ,r y ])«[x,x ; .]) ^ (2.4) 

After integration over the fermions one encounters a Pfaffian which on generic scalar field 
backgrounds is complex. This renders simulation of the model problematic. Because of this we 
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have initially focused on a related four supercharge model corresponding to the dimensional reduc- 
tion of jY = 1 super Yang Mills in four dimensions. This model has received extensive attention 
in the literature as it may be discretized on a lattice while preserving part of the supersymmetry al- 
gebra exactly 1 [||]. The J2 = 4 model has an action similar to that shown above with the restrictions 
that there are now just 3 scalars and the fermions are 4 component Majorana fields with Yukawa 
couplings related to the corresponding four dimensional Dirac matrices in Majorana representation. 

3. Lattice Theory 

The transition to a lattice theory is potentially delicate. Naive discretizations generically break 
supersymmetry at the classical level leading to the appearance of supersymmetry breaking coun- 
terterms in the quantum effective action. Some of these may be relevant which leads to a fine 
tuning problem as their associated couplings must then be fine tuned to recover supersymmetry in 
the continuum limit. 

However, super Yang-Mills quantum mechanics is, of course, super-renormalizable, which 
implies that the continuum theory contains only a finite number of superficially U.V divergent 
Feynman graphs and they occur for small numbers of loops. It is only these terms that can generate 
new supersymmetry violating interactions on the lattice. This is the content of a powerful theorem 
on lattice Feynman diagrams due to Reisz [^] but it is easy to understand - the lattice propagators 
resemble continuum propagators for light modes and only depart strongly for modes near the cut- 
off. In a theory with no cut-off sensitivity these U.V modes are irrelevant and the lattice diagram 
converges for vanishing lattice spacing to its continuum counterpart. 

In our quantum mechanics case there is only one potentially U.V sensitive diagram - a one 
loop fermion tadpole correction to the scalar field. However the gauge index structure of this term 
ensures it must vanish even in the lattice theory. Thus we conclude that any discretization of the 
continuum theory prohibiting doublers will flow to the correct supersymmetric continuum limit 
without fine tuning. 

Hence we have employed the following lattice action: 




Here ^W,^ 2 ) label the upper and lower components of the fermion field. The matrices T, are 
essentially the T, up to an redefinition of the Euclidean time direction in order that the lattice 
kinetic term takes the explicitly antisymmetric form given above 2 
The lattice covariant derivative is just 

D+f = U!(t)f(t+i)ul(t)-f(t) (3.2) 
'A property it shares with the £2 = 16 model 

2 This trick happens automatically in the twisted formulations described in [g], In the i? = 4 case described here it is 
possible to rewrite the Pfaffian as the determinant of matrix operator of half the size and discretizations of this operator 
can then choose a kinetic term proportional to the unit matrix. We have also used this equivalent formulation in our 
numerical work 
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and the lattice coupling is given by K = ^M- where we measure the inverse temperature j8 in units 
of the 't Hooft coupling A and M is the number of lattice points. The continuum limit is simply 
M — > oo at fixed /3. We also look for limiting behavior for large N. 

4. Simulation details 

It is not hard to show that the Pfaffian arising after integration over the fermions is real and 

i 

positive semi-definite allowing it to be replaced by the term (M^M) 4 where M is the fermion 
operator described in the previous section and we use antiperiodic thermal boundary conditions for 
the fermions. This latter weight is generated by an auxiliary integration over pseudofermion fields 
F,F' with action 

•Spf =F t (M t M) _3 F (4.1) 

Following the usual RHMC algorithm the latter is approximated by a partial fraction expansion of 
the form 

1 Q a 

l^+L^V (4 - 2) 

where the a,,/3; are computed using the remez algorithm [Q]. The resultant action is simulated 
using standard HMC techniques together with a multimass CG-solver to solve the Q linear systems 
in a time independent of Q [§]. 

We have simulated lattices with 5 12 points at values of the inverse temperature j3 ranging 

from 0.01 3.0. Typically we have amassed between 10 3 — 10 4 % = 1 HMC trajectories for 

systems with a number of colors in the range N = 5 16. 

Our primary observables are the absolute value of the trace of the Polyakov loop P, its as- 
sociated susceptibility %p = {<P y — < P > 2 ) and the mean energy of the system E. A simple 
scaling argument shows that the latter is simply related to the expectation value of the bosonic ac- 
tion Sb- We have also examined the mean extent of the scalar fields given by W = -V / dlx\ji.\P (jJ.). 

We have simulated both the full supersymmetric theory and a quenched theory in which the effects 
of dynamical fermions are removed. 

5. Results 

Our results for the Polyakov line P in both dynamical and quenched ensembles are shown in 
figure | We see good evidence for 't Hooft scaling in both cases. In the quenched theory we see 
evidence of a rapid crossover between a low temperature confined phase and a high temperature 
deconfined phase. This crossover seems to strengthen with increasing number of colors N. How- 
ever no such crossover is seen in the supersymmetric case and the system seems to exist in a single 
deconfined phase. This conclusion is strengthened by looking at the susceptibility %p m figure || 
which shows a a sharp peak around j3 c ~ 0.9 which grows with ,/V in the quenched theory but no 
such peak in the supersymmetric theory - a very broad peak is seen at small ,/V and for large /3 but 
this seems to move rapidly to larger /3 with increasing ,/V and we conclude it will not survive the 
large ,/V limit. Further evidence of this dramatic difference between quenched and supersymmetric 
theories is seen in the energy E shown in figure [| In the case of the quenched theory the linear 
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Figure 1: Polyakov line - quenched and supersymmetric cases 
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Figure 2: Polyakov susceptibility - quenched and supersymmetric cases 
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Figure 3: Mean energy - quenched and supersymmetric cases 



slope visible in j3 < E > for large j8 yields the non-zero vacuum energy of the model while the 
horizontal regime below j3 c is consistent with classical equipartition and the appearance of 0(N 2 ) 
weakly coupled degrees of freedom. The dynamical system is quite different - we see good 't Hooft 
scaling but only weak j3 dependence - the large /3 behavior is qualitatively similar to the quasi free 
behavior seen at high temperature. Again the simplest interpretation is that the deconfmed phase 
persists for all finite temperatures. 



6 



Gauge-gravity duality 



Simon Catterall 



6. Conclusions 

Holographic dualities offer the possibility of studying string and quantum gravity theories us- 
ing the tools of (supersymmetric) gauge theory Typically the gauge theories that arise in these 
contexts are strongly coupled which means that lattice methods are useful. Furthermore, low di- 
mensional examples exist which allow for naive discretization methods which nevertheless regain 
supersymmetry with little or no fine tuning. These models can be simulated relatively easily using 
modern algorithms even in the 't Hooft large N limit. 

In this talk we have presented results for a toy model with just £2 = 4 supercharges. Many 
features of this model are thought to be common to the target £2 = 16 model but the toy model is 
easier to simulate primarily because it possesses a real positive definite Pfaffian after integration 
over the fermions. 

We have contrasted the supersymmetric simulations with quenched simulations. In both cases 
good large N scaling is seen. In the quenched case the system exists in two phases in the limit 
N — > oo - a low temperature confined phase with non-zero vacuum energy and a high temperature 
deconfined phase with 0(N 2 ) quasi free degrees of freedom. In the supersymmetric case the de- 
confined phase appears to extend to low temperatures. This is in agreement with conjectures about 
the J2 = 16 supercharge theory [Q]. The latter theory has recently been studied using momentum 



space methods in [10] with interesting results. Perhaps most encouraging is the observation that 



the phase fluctuations of the Pfaffian were rather small over the temperature regime studied. We 



are currently investigating the J2 = 16 supercharge theory and hope to report results soon [JTT ]. 

SC is supported in part by DOE grant DE-FG02-85ER40237. TW is partly supported by a 
PPARC advanced fellowship and a Halliday award. Simulations were performed using USQCD 
resources at Fermilab. 
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